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ABSTRACT 

We present three-dimensional simulations on a new mechanism for the detonation of a sub-Chandrasekhar 
CO white dwarf in a dynamically unstable system where the secondary is either a pure He white dwarf or a 
He/CO hybrid. For dynamically unstable systems where the accretion stream directly impacts the surface of 
the primary, the final tens of orbits can have mass accretion rates that range from 10~ 5 to 10~ 3 M Q s" 1 , leading 
to the rapid accumulation of helium on the surface of the primary. After ~ 10~ 2 M Q of helium has been ac- 
creted, the ram pressure of the hot helium torus can deflect the accretion stream such that the stream no longer 
directly impacts the surface. The velocity difference between the stream and the torus produces shearing which 
seeds large-scale Kelvin-Helmholtz instabilities along the interface between the two regions. These instabilities 
eventually grow into dense knots of material that periodically strike the surface of the primary, adiabatically 
compressing the underlying helium torus. If the temperature of the compressed material is raised above a crit- 
ical temperature, the timescale for triple-a reactions becomes comparable to the dynamical timescale, leading 
to the detonation of the primary's helium envelope. This detonation drives Shockwaves into the primary which 
tend to concentrate at one or more focal points within the primary's CO core. If a relatively small amount 
of mass is raised above a critical temperature and density at these focal points, the CO core may itself be 
detonated. 

Subject headings: accretion, accretion disks — instabilities — novae, cataclysmic variables — white dwarfs 
— binaries: general — supernovae: general 



1. INTRODUCTION 

White dwarfs (WDs), the end point of stellar evolution for 
most stars, are extremely commo n, with about 10 10 of them 
residing within the Milky Way dNapiwotzkil 120091) . WDs 
are frequently observed in binary systems with normal stel- 
lar companions, and, less frequently, with compact stellar 
companions. Double degenerate (DD) systems are those in 
which the companion is another WD. Typica lly, DD systems 
are fo rmed via common envelope evolution dNelemans et alj 
12001b . Often, the final result of this evolutionary process is a 
binary consisting of a carbon- oxygen (CO) WD and a lower- 
mass helium WD companion (Napiwot zki et al.ll2.Q07l) . 

For decades, Type la supernovae (SNe) have been em- 
ployed as standard candles. Still, even the preferred mech- 
anism exhibits su bstantial variability dTimmes et al.l 120031 : 
iKasen et"aT] 120091) . Complicating this issue further is that 
WDs may not need to exceed the Chandrasekhar limit to be 
capable of exploding; other mechanisms include collisions 
between WDs in den se stellar environment s such as the cores 
of globular clusters (Ros swog et al.ll2009al). tidal encounters 
with moderately massive black holes (iRosswog et aD I2008L 



I2009bt iRamirez-Ruiz & Rosswogll2009D . or as we are intro- 
ducing in this paper, the accretion of dense material from a 
He WD companion. 

At the distance of tens of WD radii, gravitational radiation 
can bring two WDs close enough for the less-massive WD 
(secondary) to overfill its Roche lobe and begin transferring 
mass to the more-massive WD (primary). Because of the 
mass-radius relationship of WDs, mass tr ansfer is often un- 
stable and can eventually lead to a merger dMarsh et al.ll200l: 
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iGokhale et al.l 12007). For binary mass ratios close to unity 
the circularization radius drops below the radius of the pri- 
mary Z?i, and thus the accretion stream will directly impact 
the primary's surface. 

Most studies of dynamically unstable DD systems focus on 
the evolution of the post-merger object and, assuming that the 
final object has a mass larger than the Chandrasekhar limit, on 
how the merged remnant may e ventually lead to a Type la SNe 
(lLividl2000l:lYoon et al.ll2007l) . These models all assume that 
the rapid accretion that precedes coalescence is uneventful. In 
this Letter, we present three-dimensional hydrodynamics sim- 
ulations that demonstrate that explosive phenomena are an in- 
evitable results of the high accretion rates characteristic of the 
final stages of a merger in a DD system. In some cases, these 
explosive phenomena may lead to the complete detonation of 
the CO primary. 

2. NUMERICAL METHOD AND INITIAL MODELS 

The lead-up to the merger has only been studied re- 
cently due to the sensitivity of the accretion to the trans- 
fer of angular momentum. Several groups have simulated 
these mergers in th r ee dimensions using par ti cle-based codes 
dBenz et alJ 119901: iRasio & Shaoirol | 1995b ISegretain et alJ 
Il997t iGuerrero et alJl2004t iDan et aljl2008t 120091) . with the 
latter being able to resolve the final ~ 100 orbits of the binary 
prior to coalescence. SPH simulations, with their excellent 
angular momentum conservation, are well-suited to determine 
the binary orbital evolution. However, because the SPH par- 
ticles have a fixed mass and typical accretion rates are > 10~ 5 
the total mass of the system per orbit, the accretion stream is 
only resolved by hundreds of particles at best. Consequently, 
SPH approaches are unable to resolve any of the accretion 
stream's detailed structure, including how the stream impacts 
the surface of the primary. 

Eulerian codes are able to provide high resolution based 
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FIG. 1 . — Setup of run Ba in FLASH, with several annotated regions of interest. The color scheme shows log T through a slice of the orbital plane. The cyan 
circle shows the spherical outflow boundary condition centered about the secondary, while the cyan wedge shows a cross-section of the cone used as the mass 
inflow boundary. The dashed contour shows the system's Roche surface. The white boxes with labels are described in Sectionf5] 



on any combination of local and global simulation proper- 
ties, and thus are ideal for systems where the regions of 
interest do not overlap with the regions o f highest density 
dNew & Tohlind Il997t ISwestv et all 120001) . Unfortunately, 
it is difficult to accurately simulate the overall orbital evo- 
lution for many orbits in grid-based codes due to the non - 
conservation of angular momentum (iRrumholz et al.l 120041) . 
although a c areful choice of coordinate system has led to re- 
cent success dD'Souza et al.lfcOO&rMofl et al.ll2007l) . Because 
the stability of the system depends on both the structure of 
both the primary and the secondary, the grid-based simula- 
tions that investigate stability must fully resolve both stars, 
and thus the stream and the primary are not very well-resolved 
in these simulations. 

Here we follow a hybrid approach that combines the 
strengths of both SPH- and grid-based methods. We follow 
the impact of the mass transfer onto the orbital evolution by 
performing an SPH simulation (Dan et al. 2009, in prep.) and 
use the obtained, time-dependent results for orbital separa- 
tion a nd M as boundary conditions for FLASH dFryxell et al.l 
120001) . This combination allows us to realistically determine 
the dynamic stability of the binary while not sacrificing our 
spatial resolution of the accretion itself. 

In FLASH, the primary is built under the assumption of 
spherical hydrostatic equilibrium and a constant tempera- 
ture of 5 x 10 5 K, and is mapped into the simulation ex- 
plicitly. The secondary is not mapped directly into the grid 
and is taken to be a point mass surrounded by a spherical 
outflow boundary condition. The simulation is performed 
in a non-inertial frame rotating with angular frequency u> = 
\J G(M\ +M2)/a 3 about the barycenter. The primary's grav- 
ity is calculated using a multipole expansion of its potential, 
and because the quadrupolar contribution of the secondary is 
negligible, we approximate the secondary's gravitational field 
as a monopole. Because the choice of frame and the pres- 
ence of the secondary changes the hydrostatic configuration 



of the primary, we first allow the primary to slowly relax in 
the presence of the secondary's potential by removing a frac- 
tion of its kinetic energy every time-step over many dynamical 
timescales. 

The nuclear composition of the fluid is ev olved via the 1 3 
element a-chain network of iTimmesI ( 119991) . Pressures and 
temp eratures are calculated us ing the Helmholtz equation of 
state dTimmes & Swestyll2000l) . The accretion rate provided 
by the SPH calculation is used to generate a quasi-hydrostatic 
inflow boundary condition to replicate the flow of matter from 
the secondary to the primary (Figure [TJ. The boundary we 
use is cone-shaped with the apex located at LI, and mass is 
added to the cone under the assumption that the stream is in 
hydrostatic equil ibrium in the directio n perpendicular to the 
barycentric line dLubow & Shdfl976l) . The flow is directed 
towards the barycenter of the system at velocity v = c s ~ 10 8 
cm/s, the sound speed within the secondary. In all simulations 
presented here, the stream is assumed to be pure He. The four 
simulations that we have performed are summarized in Table 

m 

3. ACCRETION STREAM INSTABILITIES 

In a binary system undergoing mass transfer, the fate of the 
accretion stream is determined by ratio Ri/Rh- For all three 
of the binaries we describe R^ <Ri and thus the stream would 
directly impact the surface of the primary assuming the stream 
was collisionless. However, the stream can eventually be de- 
flected by the ram pressure of the thick torus of helium that 
gradually accumulates on the primary's surface, preventing a 
direct impact. 

As the sound speed within the torus is far smaller than the 
Keplerian rotation velocity, the collision between the torus 
and the stream is highly supersonic (M > 10), which leads to 
the development of a standing shock (FigureQ] region I). This 
shock establishes pressure equilibrium across the stream-torus 
interface. Because of the temperature increase of the material 



Surface Detonations in DD Systems Triggered by Accretion Stream Instabilities 



3 



TABLE 1 
Simulation results. 



Run 


^evol 


1 . 2 

( min 


Primary 


Secondary 


^torus,max^ 


Surf. 


Mass ejected 


Burning yield (10 3 M@) 




(s) 


(10 6 cm) 


M(M e ) 


Type 


M(M Q ) 


Type 


(10- 2 M o ) 


det. 


(1O- 2 M ) 


20 < A < 32 


A > 36 


A 


1885 


4.9 


0.67 


CO 


0.45 


He 


14 


No 








B 


434 


3.7 


0.9 


CO 


0.45 


He 


9.7 


Yes 


5.9 


11 


1.9 


Ba 4 


25 


1.8 


0.9 


CO 


0.45 


He 


11 


No 








C 


335 


3.4 


0.9 


CO 


0.6 


He/CO 


5.0 


Yes 


6.0 


16 


5.6 



"Total evolution time. 
^Minimum grid cell size. 

^Maximum mass of the helium torus accumulated during the run; this is either 
rf Run Ba is initialized from a checkpoint produced by ran B at t = 403s. 

falling from LI to the surface of the primary, the torus consists 
of material heated to a temperature close to the virial temper- 
ature of GMim p /R\k\, ~ 10 8 K. Radiation pressure becomes 
competitive with degeneracy pressure at these temperatures 
and densities, and as a result the torus has a density smaller 
than the incoming accretion stream. 

Because the accretion stream and the helium torus do not 
move in the same direction, the difference in velocity at the in- 
terface is substantial, U\-Uo^ \J GM\ jr\ (y/2— cos 6), where 
9 is the angle between the stream and torus velocities. This 
shear can lead to the rapid growth of Kelvin-Helmholtz insta- 
bilities along the interface between the two regions (Figure[T| 
region II). An important factor in determining which pertur- 
bations are unstable is the magnitude of the restoring force 
perpendicular to this interface. The a cceleration g provides 
stability for perturbations with k < (IChandrasekharl 1 96ll) 



the value just before detonation or at the end of the ran if there was no detonation. 
0.7 



_ g(ai-a 2 ) 
a\a2(U\ -U2) 



(1) 



where a,\ = p\/(p\+pi) and = pi/(p\+ pi)- As long as 
g > 0, perturbations with wavelength A > A max = 27r/£ m j n will 
not grow. Within the primary's Roche lobe, g is determined 
by the combination of the primary's gravity and the Coriolis 
force 

oom GMi\riXv\ 
g = 2n\v]-— - ^ j ^ . (2) 



in 



1*1 M 



Because the Coriolis force is always perpendicular to the ve- 
locity, its contribution to g is always equal to the magnitude 
of the force. Only the component of the primary's gravity that 
is perpendicular to v*will contribute to g. Near LI where the 
cumulative gravitational forces are parallel to the velocity, the 
Coriolis term is dominant and is directed away from the pri- 
mary, and thus g is positive. As the stream falls toward the 
surface of the primary, U\—U2 continually increases, and the 
restriction on k is progressively less limiting. And because 
the Coriolis effect is oc | v| oc \r\ \~ l l 2 while the primary's grav- 
ity is oc |Fi|" 2 , the primary's gravity eventually becomes the 
dominant contributor to the restoring force. As the primary's 
gravity is opposite to the Coriolis force, g eventually changes 
sign, at which point all wave modes become unstable (Figure 
Q] region III). 

The helium torus, which is both clumpy and has a radial 
density profile, drives the seed vertical displacements in the 
accretion stream. The torus is highly turbulent, with the 
largest bulk motions being driven at scales comparable to 
the width of the accretion stream (<~ 10 8 cm) and at speeds 
that are highly supersonic. Therefore, the power spectrum 
of the resulting turbulent cascade is given by Burgers turbu- 
lence (E oc k~ 2 ). Consequently, perturbations with large wave- 




FlG. 2. — The hashed region on the lower right shows systems that are def- 
initely dynamically stable, while the sticked region in the upper left shows 
systems that are definitely unstable. The thick solid line shows the transition 
from direct impact to disk accretion, while the dot-dashed curve shows the 
separ ation between sub- and super-Eddington accretion in the weak coupling 
limit (Marsh et al. 2004). The contours show where T<i yn = T3 Q for differ- 
ent stream densities p s 5 = p s /10 5 assuming a torus density of 3 X 10 5 g 
citT 3 . Surface detonations are likely for systems that accrete with p s for 
which Td yn < T3 Q and then experience an increase in p s until rj yn > T} a . 

lengths are seeded with large initial amplitudes. Addition- 
ally, power is transferred from large k to small k by the tidal 
stretching of waves as they fall in the primary's gravity. A nu- 
merical calculation of the growth rate of the instabilities show 
that even the largest unstable modes, which have the slowest 
growth rates, experience a ten-fold increase in amplitude in 
the time it takes to fall from LI to the primary's surface. Be- 
cause the flow is three-dimensional, instabilities larger than 
the stream size cannot form, and the end result is that the 
stream consists of knots of material that are comparable in 
size to the stream itself (FigureQ]). 

4. SURFACE DETONATIONS 

We find that instabilities in the accretion stream are present 
in all four runs once the helium torus has acquired enough 
mass to deflect the stream. The material that is compressed 
against the surface of the primary by the dense knots formed 
in the stream can be heated to billions of degrees Kelvin (Fig- 
ure [U region IV). At the onset of mass transfer, the density 
of the stream is too small to heat the underlying material ap- 
preciably. This allows ~ 1- 15% of a solar mass of helium to 
accumulate on the surface of the primary. As the density of 
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FIG. 3. — Volumetric snapshots showing the temperature T and logpne = pXne during the explosive event on the surface of the primary in run C. A dense knot 
of helium that formed through Kelvin-Helmholtz instabilities in the accretion stream is marked by a magenta arrow. The evolution of the surface detonation is 
described in Section[3] 



the accretion stream increases with the decreasing separation 
between the WDs, the temperature of the compressed layer 
can exceed the temperature required to explosively ignite he- 
lium via the triple-alpha process. This leads to a gravitation- 
ally confined surface detonation that propagates through the 
helium torus that wraps around the surface of the primary at 
approximately the Keplerian velocity (Figure O. 

The conditions required for a successful surface detonations 
are summarized in Figure [2] Only systems that are dynami- 
cally unstable will ever have accretion rates capable of ignit- 
ing surface detonations via stream instabilities, and thus the 
lower-right region of the figure is excluded. The question of 
whether the surface detonation will occur in a given binary 
is tricky to answer because the process is highly stochastic; 
the impact trajectory, density, and geometry of the knots that 
initiate the detonations vary considerably over the course of 
the simulations. Nevertheless, we can estimate when surface 
detonations will occur by comparing the tim escale for triple- 
q reactions T3 Q dCaughlan & Fowlerll 19881) to the dynamical 
timescale Td yn ~R\/ v eS c,i at the surface of the primary. 

For a surface detonation to involve a significant amount of 
mass, the system must accrete material at a rate for which 
i~3a > Tdyn, and then transition to a rate where T3 a < Td yn . 
By solving the equations of motion for a test particle re- 
leased from Li, we estimate the two timescales in Figure|2]for 
various stream densities for different primary and secondary 
masses. The conditions at the He/CO interface are given by 
dP/dr= d<f)/dr, where <f> ~ GM\jr\ -rif^ or and Sl toI is the 
angular velocity of the torus, which is related to the compo- 
nent of the accretion stream velocity parallel to the primary's 
surface. Conversely, the ram pressure P ram oc pv 2 applied by 
the accretion stream is dependent on the component of veloc- 
ity perpendicular to the primary's surface. 

The density structure of the torus is approximately fixed as 
the system evolves towards merger because (f> only increases 
slightly due to mass accretion. Meanwhile the stream density 
continually increases, resulting in the contour of T3 Q = Td yn 



moving leftwards in Figure[2] When T3 Q < Td yn , the increase 
in temperature of the torus can no longer be controlled by 
adiabatic expansion, which results in thermonuclear runaway 
and a detonation in the helium torus. 

5. RESULTS AND DISCUSSION 

Surface detonations are present in both of the longer runs 
with a 0.9 M Q CO primary, although run B detonates much 
later in its evolution than run C. The geometry and evolution 
of the detonations are similar in both runs. A small region 
of the surface helium is heated to a temperature of ~ 2 x 10 9 
K, leading to a detonation front that expands outwards from 
the ignition site along the primary's surface (Figure [3] sec- 
ond column). Because the helium layer is toroidal rather than 
spherical, the front runs out of fuel as it propagates towards 
the primary's poles. Consequently, the detonation splits into 
two fronts that run clockwise and counterclockwise around 
the equator of the primary (Figure [3] third column). These 
detonation fronts eventually run into each other along a longi- 
tudinal line that is opposite to the original ignition site (Figure 
[3] fourth column). The highest temperatures are produced in 
this convergence region. 

For run A, the primary's surface gravity is too low to com- 
press the helium layer above the critical temperature neces- 
sary for thermonuclear runway. No detonation was observed 
in run Ba, despite being initialized using a checkpoint from 
run B just prior to the observed detonation in B. This is be- 
cause of the stochastic nature of the ignition mechanism — 
the particular dense knot of material that led to the surface 
detonation in B has a different shape in Ba, and did not have 
a favorable geometry for igniting the helium torus. And be- 
cause Ba has twice the linear resolution of run B, we could 
only afford to evolve it for a short period of time. 

If a helium surface detonation occurs as the result of stream 
instabilities in a DD system, the resulting t ransient event 
could potential resemble a dim Type la SNe dBildsten et alj 
120071: IPerets et alj 120091) . The results of our simulations 
shows a large degree of variability which ultimately depends 
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FIG. 4. — Maximum temperature T in a mass shell as a function of inte- 
rior mass M(< r) and time ; for run B. Each shell lies on a surface of equal 
gravitational potential. The shading is logarithmically binned, with red cor- 
responding to and blue corresponding to T m j n = T max /10' 5 . The dashed 
line shows the region where Xn e > 0.02, while the dotted line shows where 
X Nc > 10" 3 . For convenience, we include a zoomed view of the surface det- 
onation as an inset. 

on the geometry of the dense knots when they strike the sur- 
face of the primary in the impact zone. However, if the post- 
detonation temperature is not much larger than 2 x 10 9 K or 
the geometry of the detonation fronts are not favorable, it 
is possible to only synthesize intermediate-mass a-elements 
(Table[Q. 

The prospect of igniting the CO core itself is attractive 
because the typical densities of moderate-mass WDs (M ~ 
1.0M Q ) are ~ 10 6 - 10 7 g cm" 3 , com parable to the de nsities 
found in t h e DD T for Type la SNe (iKhokhlov etai1ll997l) . 
iFink et all (120071) (hereafter FHR) show that when 0.1M Q 
of helium is ignited at the He/CO interface of a moderate- 
mass WD, strong shocks are launched deep into the CO core, 
which converge at a focusing point. This focusing produces 
a region where the temperature and density meet the criteria 
for explosively igniting carbon (iNiemeyer & Woosley! 1 19971 



iRopke et alJ2007l) . thus leading to a full detonation of the CO 
core. 

The linear resolution of our simulations is substantially 
coarser than the highest-resolution two-dimensional simula- 
tions of FHR, and we are certainly under-resolving the true 
density and temperature peaks. However, FHR's results ap- 
pear to be rather optimistic when applied to our model for the 
following reasons. First, FHR assumes perfect mirror sym- 
metry, which dramatically increases the amount of focusing 
by forcing the shocks to converge at a single point rather than 
a locus of points. Second, FHR does not use a nuclear net- 
work and makes the simplifying assumption that the entire 
He envelope is burned to Ni, which naturally overestimates 
the energy injected into the CO core when compared to the 
more realistic case where burning is incomplete. 

Because of these differences, the conditions for double det- 
onation are still not met even in our run with the most favor- 
able shock geometry and strongest surface detonation (run C). 
In the focusing region of run C, the temperature and density 
are 2.5 x 10 s K and 2.0 x 10 7 g cm" 3 , respectively. Conditions 
for CO core detonation may be more propitious in systems 
containing a slightly more massive primary, but ultimately the 
successful detonation of the core depends on the precise dis- 
tribution of temperature and den sity within the high-pre ssure 
regions of the convergence zone (ISeitenzahl et al.ll2009l) . 

The mechanism reported here for the detonation of a sub- 
Chandrasekhar CO WD in a dynamically unstable binary is 
not tied to a particular mass scale and therefore allows for 
considerably more diversity. As outlined above, mass transfer 
between a pure He WD or a He/CO hybrid and a CO WD be- 
fore the merger provides a novel pathway to ignite CO WDs. 
Even if a critical amount of mass is not raised above the con- 
ditions required for CO detonation, a peculiar underluminous 
optical transient should signal the last few orbits of a merging 
system. 
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